Generalized drift-diffusion model for miniband superlattices 
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A drift-diffusion model of miniband transport in strongly coupled superlattices is derived from 
the single-miniband Boltzmann-Poisson transport equation with a BGK (Bhatnagar-Gross-Krook) 
collision term. We use a consistent Chapman-Enskog method to analyze the hyperbolic limit, at 
which collision and electric field terms dominate the other terms in the Boltzmann equation. The 
reduced equation is of the drift-diffusion type, but it includes additional terms, and diffusion and 
drift do not obey the Einstein relation except in the limit of high temperatures. 
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In recent years, nonlinear charge transport in semi- 
conductor superlattices (SLs) has blossomed as a field, 
driven by the availability of many experimental results 
and by theoretical analyses and simulations of rate equa- 
tion models^^. As it often happens, these models have 
not been derived from more fundamental "first princi- 
ples formulations" such as kinetic theory. The situation 
is different depending on whether the SLs are weakly or 
strongly coupled. In weakly coupled SLs, neighboring 
quantum wells are separated by "thick" barriers and ver- 
tical transport occurs via sequential resonant tunneling 
through them. Provided intersubband scattering is much 
faster than escape times from a quantum well and the 
latter are much smaller than dielectric relaxation times, 
electrons are at local equilibrium in the subband of low- 
est energy^. Then the tunneling current density across 
a barrier under stationary conditions can be calculated 
from "first principles" using the Transfer Hamiltonian 
method^, Green functions for a SL under a constant ex- 
ternal field 3 , etc. This tunneling current (that depends 
on the electron density in the two quantum wells sepa- 
rated by the barrier and on the local value of the elec- 
tric field) is then inserted in a discrete rate equation 
model including charge continuity and a discrete Poisson 
equation^. No derivation of this very reasonable model 
seems to be known to this date, although its validity has 
been corroborated by numerous experiments. 

In strongly coupled SLs, barriers are "thin" , minibands 
are wide, and quantum wells cannot be considered as 
separate entities. Practical models to analyze nonlin- 
ear transport are of the drift-diffusion type^ or hydrody- 
namic models&t Drift-diffusion equation (DDE) models 
typically use a drift velocity obtained from a simplified 
kinetic equation and a diffusion coefficient that obeys the 
Einstein relation^. The resulting model is a variant of the 
well-known Kroemer DDE for the Gunn effect in bulk 
n-GaAA For the large fields involved and for the non- 
parabolic SL miniband energy, using an Einstein relation 
to figure out the diffusion coefficient is questionable^ 



and, in fact, incorrect except in a particular limit. Hy- 
drodynamic models are considerably more complicated 
and have been solved numerically, but not many analy- 
ses of them have been carried out. The same applies to 
quantum diffusion theories^. "First-principles deriva- 
tions" typically solve a kinetic equation numerically or 
approximately assuming a constant applied electric field 
and ignoring space and time dependence. Then a current 
density across the SL is calculated for different values 
of the field, and a drift velocity and a diffusion coeffi- 
cient are figured out. These functions are then inserted 
in a DDE. Results that are "valid for any type of SL" 
typically mean that stationary, space-independent solu- 
tions of a sufficiently general kinetic equation have been 
found numerically 3 . Again the crucial derivation of a 
rate equation model from a kinetic equation is missing. 
In this paper, we provide such a derivation starting from 
a simple Boltzmann-Poisson system that describes one- 
dimensional (ID) electron transport in the lowest mini- 
band of a strongly coupled SL: 
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Here I, e, f, n, Nrj, fcs, T, F, m* and e > are the 
SL period, the dielectric constant, the one-particle dis- 
tribution function, the 2D electron density, the 2D dop- 
ing density, the Boltzmann constant, the lattice tem- 
perature, minus the electric field, the effective mass of 
the electron, and minus the electron charge, respectively. 
The first collision term represents energy relaxation to- 
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FIG. 1: Electron density versus chemical potential at 10 K 
(M = 25.59/fcsT) and at 300 K (M = -0.45/tsT). 



wards a ID effective Fermi-Dirac distribution f FD (k;n) 
(local equilibrium)^ with collision frequency v e . The 
second collision term accounts for impurity elastic colli- 
sions: f^M^^k'WSity-Sik'MfW-fmdk' = 
24, (x,k,-k)[f(-k) - f(k)]/(Alsmkl) = Vi[f{-k) - 
f(k)]/2, provided we use the tight-binding miniband dis- 
persion relation, £(k) = (A/2) (1 — cosfcZ) (A is the mini- 
band width), and ignore transversal degrees of freedom 11 . 
For simplicity, v e and Vi will be fixed constants. 

Exact and Fermi-Dirac distribution functions have the 
same electron density, thereby preserving charge conti- 
nuity as in the BGK (Bhatnagar-Gross-Krook) models 
of collision processes^. Then the chemical potential /i 
depends on n and is found by inverting the exact rela- 
tion (PJ; cf. Fig. ^ BGK collision terms with a Boltz- 
mann distribution function, the Boltzmann limit of Eq. 
(@J, were introduced by Ignatov et al. who adapted 
collision models by Ktitorov et al. containing inelastic 
energy relaxation and elastic impurity momentum relax- 
ation termsii 

To derive a reduced balance equation for n, we shall 
assume that the electric field contribution in Eq. is 
comparable to the collision terms and that these terms 
dominate the other two. This is the so-called hyper- 
bolic limit, in which the ratio of df jdt or v{k)d f / dx to 
(eF/h) df/dk is of order e 1. Let v M and F M be elec- 
tron velocity and field scales typical of the macroscopic 
phenomena described by the sought balance equation; for 
example, let them be the positive values at which the (ze- 
roth order) drift velocity reaches its maximum. In the hy- 
perbolic limit, the time to it takes an electron with speed 
vm to traverse a distance xq = sFmI/ (e-/Vo), over which 
the field variation is of order Fm, is much longer than the 
mean free time between collisions, v~ l ~ h/(eFMl) = t\. 
We therefore define e = t\jt a — TlvmNd/ (sF M l 2 ) and 
formally multiply the two first terms on the left side of 
© by e J£ After obt aining the number of desired terms, 
we set e = 1. The solution of Eq. for e = is 
straightforwardly calculated in terms of its Fourier co- 



efficients as /(°)(fc;n) = /] °'e^ kl , with /j 
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F M = h\fv e {ve + Vi)/{el), and r e = \J{v e + Vi)/v e . 
Since f FD is an even function of k, its Fourier coefficient 
f FD is real. Note that Eq. © implies f^ 0) = f FD = n. 

We shall derive a reduced balance equation for the elec- 
tron density by using the Chapman-Enskog ansatai^: 
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The coefficients f^ m \k; n) depend on the 'slow variables' 
x and t only through their dependence on the electron 
density and the electric field (which is itself a functional 
of n). The electron density obeys a reduced evolution 
equation JfJJl in which the functionals N^ m \n) are cho- 
sen so that the / (m) (fc; n) are bounded and 27r/Z-periodic 

in k. Moreover the condition, l_^n f^ m \k;n) dk = 

2irf( m) /l = 0, m > 1, ensures that f( m \ m > 1, 
do not contain contributions proportional to the zero- 
order term /(°). N^{n) can be found by integrating 
^ over k, using ©, and inserting JSJ in the result: 

ATM(n) = -l(d/dx)J* / J, l v(k)f( m '>dk/(2ic). Then, in- 
tegration of @ over x yields 
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v(k) f i - m \k;n)dk = J{t), (7) 



where J(t) is the total current density. To find the equa- 
tions for f( m \ we insert (j^J and 10 in and equate 
like powers of e: 
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and so on. We have defined Cu(k) = eFh~ 1 du(k)/dk + 
[v e + i>i/2)u(k) + Viu{—k)/2, and the subscripts and 1 
mean that dn/dt is replaced by N^(n) and by N^^n), 
respectively. 

The linear equation Cu — S has a bounded 2ti/1- 
periodic solution provided S dk = 0. This solvabil- 

ity condition together with Eqs. ©, etc. also yield 
the previously found and the reduced equation (7J| . 

Keeping only the leading order terms in 0, we obtain 
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provided M{n/N D ) = I 1 (f ll )I {M)/[l 1 {M)l (n)], A = 
fi/{k B T), and 5 = A/(2k B T). Using Eq. ©, the di- 
mensionless chemical potential fx = jj,(n/N B ) is calcu- 
lated graphically in Fig. 2] as a function of n/Nrj, with 
(1(1) = M. Then we have M(l) = 1. In the Boltz- 
mann limit, M. = 1 for any n, and the electron cur- 
rent density in l|l(J|) has the usual drift form. Thus M. is 
a low-temperature, density-dependent correction to the 
usual drift current density. The drift velocity, vmV(!F), 
has the Esaki-Tsu form with a maximum that becomes 
vm ~ Alli(5) y/U^ / [4TiIo(8) y/ v e + vi\ in the Boltzmann 
limilii^ (/„ (6) is the modified Besscl function of the nth 
order). 

The first-order correction in (JJJ) is found by first solving 
|jHJ. After straightforward but lengthy calculations and 
setting e = 1, we obtain (here g' means dg/dn): 
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Here the density-dependent function M^u/Nd) = 
l2(fh)Io(M)/^Lo(fb)I\(M)] becomes simply the con- 
stant li (S)jl\ (5) in the Boltzmann limit. Despite its 
formidable appearance, the generalized drift- diffusion 
equation (GDDE) (|13fl is (in dimensionless units) a small 
perturbation of the drift equation (|l(Jfl , analyzed in stud- 
ies of the Gunn effect a long time agoi£*i£. Table [i] 
shows that the solution of the GDDE and (0 yield self- 
oscillations of the current with frequencies that agree 
with those measured by Schomburg et a&^. 

An often used DDE consists of Eq. (JTDJ (with M = 1) 
plus a diffusion term obeying the Einstein relation 8 : 
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The difference between the predictions of (|13|l and 119(1 
can be remarkable if the dimensionless parameter e is 
relatively large and the dimensionless coefficient 5 = 
A/(2k B T) is not small, as illustrated in Fig. [21 The pa- 
rameter values in this figure correspond to the 5.13-nm 
GaAs/0.87-nm AlAs SL of Ref. H for which e = 0.34: 



TABLE I: Numerical values of the oscillation frequencies 
t-'num, compared with the experimental value f cxp for five of 
the SLs of Ref. Ugl together with the corresponding applied 
voltage $. dw and ds are well and barrier widths, respec- 
tively. 
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FIG. 2: (a) Current (Jo = evMNo /I) vs. time during self- 
oscillations for a 100-period SL at 300K, as described by the 
GDDE in the Boltzmann limit (solid line) and by the DDE 
(dashed line), (b) Comparison between the dipole wave for 
the DDE, (1), and the dipole wave for the GDDE, (2). 



N D = 0.84 x 10 n cm- 2 , A = 43meV, v e = v t = 10 13 Hz, 
xq/1 = 0.75. We have selected bias and boundary con- 
ditions so that dipole mediated current self-oscillations 
occur in this SL: voltage bias divided by SL length equals 
1.2 F M , and T = 2Jl/(eN D v M ) at both SL ends. The 
difference in oscillation frequency and wave shape can be 
explained by taking into account the equal-area rule as 
in the theory of the Gunn effect: the taller wave of the 
GDDE moves at a slower average speed than the wave of 

m- 

For SLs with a smaller value of e, the difference be- 
tween the predictions of the GDDE and the DDE lj*H3|) is 
smaller. Is there a limit in which these equations agree? 
To explore this, we calculate the deviation of drift ve- 
locity and diffusion coefficient in the GDDE from the 
Einstein relation (setting n — No and fi = M): 
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FIG. 3: Ratio R(F) at: (a) 10 K, and (b) 300 K. (c) Relative 
error of the Boltzmann limit result with respect to using the 
Fermi-Dirac distribution. 
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In the Boltzmann limit, exp(/i — <5 + Scoskl) <C 1, and 
we can substitute the modified Bessel functions I S (S) 
instead of 1 S (M) in the previous formula. Moreover, 
Eqs. © and Q give n w e^ 1 ' 8 I Q (S) p N D , where po = 
m*k B T/{-Kh 2 N D ). For n = A^, the Boltzmann limit 
holds provided p ^ 1- If we also have S — A/ (2/csT) <C 
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Ignoring correcting terms, in this limit the right hand side 
of Eq. (|21|l becomes 1 and the Einstein relation holds. 
Fig. shows the deviation from the Einstein relation at 
different temperatures, either using the Fermi-Dirac dis- 
tribution in Eq. |JT]|. using its Boltzmann limit, or the 
two-term approximation (|21Jl . Deviations are more ap- 
preciable at low temperatures. In the limit k B T 3>max 
(A, Trh 2 N D /m*), the GDDE O becomes the DDE O 
up to terms of order eS if we set A = 1 + 0(e) ~ 1. 

In conclusion, we have derived a generalized drift- 
diffusion model for charge transport in miniband su- 
perlattices by means of a consistent Chapman-Enskog 
method. At all temperatures, its predictions deviate ap- 
preciably from those of the usual DDE with the Esaki- 
Tsu drift velocity and diffusion obeying the Einstein re- 
lation. The DDE holds in the limit e < 1, k B T >max 
(A,Trh 2 ND/m*), which is not very realistic for many 
strongly coupled SLs, even at room temperature. De- 
tailed analyses and comparison between the predictions 
of the two DD models and those of the original kinetic 
equation will be considered in future works. 
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